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The Beresinskii-Kosterlitz-Thouless (BKT) physics of vortices in two-dimensional superconduc- 
tors at finite magnetic field is investigated by means of a field-theoretical approach based on the 
sine-Gordon model. This description leads to a straightforward definition of the field-induced mag- 
netization and shows that the persistence of non-linear effects at low fields above the transition is a 
typical signature of the fast divergence of the correlation length within the BKT theory. 
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The Beresinskii-Kosterlitz-Thouless (BKT) transition 
[1], namely the possibility to have a phase transition with 
a vanishing order parameter but algebraic decay of the 
correlations, is undoubtedly one of the most fascinating 
aspects of collective phenomena. It finds experimental 
realizations in a wide range of systems, as superfluids or 
superconducting (SC) films[2-5] and recently cold atomic 
systems [6]. One of the key ingredients of the BKT tran- 
sition is the existence of vortices, that unbind in the high- 
temperature phase leading to an exponential decay of the 
correlations. In order to treat such unbinding transition 
a very fruitful analogy was to represent the vortices as 
charges performing a Debye-Huckel screening transition 
in a neutral Coulomb-gas problem, for which a renormal- 
ization group procedure can be implemented [2]. 

One specially interesting extension of the BKT transi- 
tion is when a magnetic field is present, which will im- 
pose a population of vortices with a given vorticity in the 
system. This has found recent experimental application 
to thin films[4, 5] or layered high-Tc superconductors[7]. 
Even in cold atomic systems, a magnetic field can be 
mimicked by imposing a rotation on the condensate [6, 8]. 
For all these systems it is thus crucial to predict theoret- 
ically how the magnetic field will affect the BKT transi- 
tion and the various physical observables. 

Due to the strong interest of such a question, this prob- 
lem has been addressed in the past [2, 9-11]. Unfortu- 
nately, contrarily to the case of the B = transition, 
the efforts have been partly unsatisfactory. In particular 
most of the literature on the subject rested on extend- 
ing the mapping to the Coulomb-gas problem, where the 
effects of the magnetic field can be incorporated as an 
excess of positive charges. However this mapping gives 
the physical observables as a function of the magnetic in- 
duction B instead of the magnetic field H, which is not 
convenient to describe the physics at low applied field. 

An alternative approach to the BKT transition, which 
is of course well known for B = 0, is to use the map- 
ping onto the sine-Gordon problem[2, 12], which was 



reviewed recently both in the context of quasi-2D su- 
perconductors [13, 14] and cold atomic systems [8]. In 
this Letter we show that this description provides a very 
simple and physically transparent way to deal with the 
finite magnetic field case. In our scheme the physical 
observables have a straightforward definition, and the 
role of both B and H is clarified. In addition we also 
present a variational calculation of the field-induced dia- 
magnetism in thin films. It leads to a detailed description 
of the Meissner phase below Tbkt and of the appearance 
above Tbkt of a non-linear magnetization at relatively 
low fields, in contrast to what expected from standard 
Ginzburg-Landau (GL) SC fluctuations [15]. 

As a starting model we consider the XY model for the 
phase of a 2D superconductor[2] 
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Here 6*^^ is the SC phase on two nearest-neighbor sites 
(j, j) of a coarse-grained 2D lattice, J — $Q(i/167r'^A^ is 
the 2D superfluid stiffness for a film of thickness d and in- 
plane penetration depth A, and we employed a minimal- 
coupling scheme for the vector potential A, with Fij = 
(27r/$o) // A-dl, and $0 = hcl2e the fiux quantum. Due 
to the periodicity of H when Q Q ^-2-^, beyond long- 
wavelength phase excitations where Qi — 9j « aV0 varies 
smoothly on the lattice scale a, vortex configurations are 
allowed where = ±27r over a closed loop. They 

emerge clearly by performing the standard dual mapping 
of the model (1)[16]. This allows us to write the partition 
function of the system as a functional integral over a 
scalar field 4> as, Z — J Vcj)! 
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where 4> depends on the in-plane coordinates r only while 
A depends in general also on the z coordinate. The 
5{z) function gives the proper boundary conditions for 
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a truly 2D case (whore there is no SC current outside 
the plane). In the physical case of a SC film of thick- 
ness d we assume that the sample quantities are aver- 
aged over \z\ < d/2. In Eq. (2) we defined K = nJ/ksT 
and g = 2ne~^'^, where /x is the chemical potential 
of the vortices and e~^^ their fugacity (/? = l/fcsT). 
Whfie in the XY model ///J is fixed, ^xy — 7r^J/2, 
we consider it as an independent variable[13]. In the 
dual representation (2) of the XY model (1) the cosine 
term accounts for vortex excitations: indeed, since (j) is 
the dual field of 9, a vortex, which is a ±27r kink in 
the 9 variable, is generated by the operator e~^^e^^'^'^ . 
At high T (j) localizes in a minimun of the cosine and 
its conjugate field 9 is completely disordered, i.e. the 
system looses the superfiuid behavior. The interaction 
V^(r) between vortices (or charges in the Coulomb-gas 
analogy[2, 12]) follows from the Gaussian part of the ac- 
tion (2), V{r) = /rf2ke*"-y(k) where l/(k) = (|<;i(k)|2), 
and it is logarithmic since V{k) = 2nK/k.'^. 

The physical observables can be easily read out from 
Eq. (2) and the free energy F = — fc^Tln Z. For example, 
the electric current is: 



action reduces to: 
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and it is purely transverse, as expected for vortex exci- 
tations. The magnetization M = (B — H)/47r is defined, 
as usual, as the functional derivative of F with respect to 
B(r, z) = Vx A. By integrating by part, we can write the 
last term of Eq. (2) as (2i/$o) / drdz [B(r, z)-i(/>(r)](5(z), 
so that: 



M(r) 



dz-— r = -z " U(r)) , (4) 
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which leads to = c(V x M)[17]. Finally, by exploiting 
the fact that e~''^e^'^'* is the operator which creates up 
and down vortices with density n± respectively, we have 
a straightforward definition of the average vortex number 

np — a?{{nj^) + {nJ)) and of the excess vortex number 
n = a'^{{n+) — (n_)) per unit cell as a function of as: 



np = 2e-'^''{cos{2(i))), n = 2e-^''(sin(2.?i)). (5) 

In Eq. (4) the average value of (f) is computed with 
the action (2), so that it gives M as a function of the 
magnetic induction B. To obtain M as a function of 
the applied field H we must use the Gibbs free energy 
g = -ksTlnZ, where Z = j VcpVAe'^ and: 
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where 1/A = d/2X^ = Sn'^KkBT/'i'l Here is 
the magnetic field generated by Jext in the vacuum, 
i.e. it satisfies the same Maxwell equation as H, but 
it is not constrained to the boundary condition that 
B = in the SC film. Thus, using the Laplace for- 
mula HO(r) = (1/c) / (iV[Je,t(r') x (r - r')]/|r - rf. 
The effect of integrating out the B field is twofold. First, 
one introduces an effective screening of the vortex po- 
tential V(r). Indeed, thanks to the kA~^ term in Eq. 
(6), V(r) ~ log(r/a) up to a scale of order A, and then 
decays as A/r[2, 10]. Second, one couples directly the 
dual field to the reference field H° used in the exper- 
iments. We thus expect that in the Meissner phase M 
includes automatically the demagnetization effects, i.e. 
— 47rM = H°/(l — ry)[17], where rj is the demagnetization 
constant which depends only on the sample geometry and 
is near to 1 in a film, 77 ~ 1 — d/R[17, 18], where R is the 
transverse film dimension. 

The model (6) and the constitutive equations (3)-(5) 
establish a clear and general theoretical framework to ad- 
dress the physics of 2D SC films in a magnetic field. To 
illustrate their usefulness we solve them by using a vari- 
ational approximation. The idea is to replace the cosine 
interaction in Eq. (6) with a mass term A^0^, where A 
is determined self-consistently by minimizing the varia- 
tional energy Qvar ~ Go + T{S — Sq), Sq being the trial 
action. At H° = a finite A appears above Tbkt, which 
signals the localization of in a minimum of the cosine, 
and cut-off at a scale 1/A the logarithmic vortex poten- 
tial V^(r). This allows for the proliferation of free- vortex 
excitations. We then consider the case of a perpendicular 
field H = Hz (in the following we drop the superscript 
0) slowly varying over the film. To account for it we 
introduce in the trial action an additional variational pa- 
rameter H , coupled linearly to / drd> in analogy with Eq. 
(6), so that only the 4>{k) component at the minimum k 
value kmin — 1/i? couples to H. A finite system size R 
is needed to have finite demagnetization in the Meissner 
phase, but its role at large fields (and in general above 
Tbkt) is negligible. The trial action is: 

^0 = ^ E [G-^(k)<^(k)0(-k)] + |iif0(k„,„), (7) 
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H satisfies the Maxwell equation V x H = {A'K/c)3ext 
for a given distribution Jf j./ of external currents. By 
integrating out A in the radial gauge V • A = 0, the 



where R^ is the film area and G~^(k) = {k'^^kh.-'^ + 
{A/a}'^)/'KK. According to Eq. (4) the magnetization is 
related to A and H as 
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kBT AttK Ha'^ 
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where M is the dimensionless magnetization and 
(AR/a)^ = + 1/RA is the intrinsic (i.e. T and H 

independent) cut-off. By minimizing Gyar with respect 
to (A, H) we derive the two self-consistent equations: 

4:Kg{A + Aa)^ cosh(M) = A^ (9) 
A^ tanh(A7) ^ nH{4:7:K) - MA% (10) 

where Aa = Aj^R/a and tih = Ho^/^q is the flux per 
unit cell. Finally, Eq. (5) leads to: 



np 



A^/AttK, n = nH - MA^/AttK. 
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We note that using Eq. (11) the two Eqs. (9)-(10) can 
be related to similar expressions derived in Ref. [9, 10]. 
Nonetheless, a clear connection to the magnetization and 
to the role of H vs B was lacking in these papers. As a 
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FIG. 1: (Color online) (a)-(b) M{H) and nF^H) below Tbkt 
from the numerical solution of Eqs. (9)-(10). The arrows in- 
dicate Hc\ according to Eq. (14). At T = QATbkt M{H) 
shows a sharp kink at Hci, where np drops abruptly to small 
values. At higher T these features are partly smoothened 
out by thermal smearing, (c) Temperature dependence of the 
mass term Ah=o. Notice that Ah=o ~ Aa at Tbkt, where 
A = 6.9 X 10^ A, Afl ~ 4 X 10"'' and Aa = 1.9 x 10"^ 



prototype of 2D system we consider a single layer of un- 
derdoped Bi2212, with J(T) = Jo(l - T/Tmf) to mimic 
the bare T dependence due to quasiparticles, Jq — 180 
K and Tmf = 120 K, which gives Tb^t = 84 K. For 
d = 15 A as the typical interlayer distance the magneti- 
zation (8) is given in units of ksT/d^o = (4.4 x 10"^T) 
G (or (4.4 X T) A/m in the notation of Ref. [7]). More- 
over, we use R ^ lO^a as the typical sample size, 
with a = 40 A, and choose fi — l.bfixY- With this 
choice of parameters one has always i? ^ A , so that 
Aa ^ 1/A > Ar - 1/VRA up to Tbkt- 

At H ^ Eq.s (9)-(10) are satisfied for Af = and 
A solution of the equation AKg{AH=o + Aa)^ — A|^^q. 
At Aa = the solution Ajj^^ = {iKg)^^^"^-^^ is finite 
only at K < 2, which identifies the BKT transition at 



K = 2 (i.e. Tbkt = t^J{Tbkt)/'2). When a finite cut- 
off Aa is introduced Ah=o approaches A a at Tbkt, and 
vanishes as Ah=o = V^T{gA^^^ as T ^ 0, giving A <C 
A\,Ar already at T < 0.9Tbkt, see Fig. Ic. Observe 
that at H = 0, where the same number of ± vortices are 
thermally- induced, n = 0, and one can parametrize np 
in Eq. (11) via the vortex correlation length ^ as 1/^^ = 
np/a"^ = A^^^/iTTKa'^. 

At H 7^ a finite M appears, which modifies also 
the A value. In general, at low field A keeps the zero- 
field value A{H) « Aj^^o and M grows linearly with 
H. By further increasing H, A grows with respect to 
Ah=o and M enters a non-linear regime. The slope of 
M vs H, the absolute value of M and the crossover field 
differ substantially above and below Tbkt- Let us first 
analyze the case T < Tbkt, i-e. K > 2. For small M one 
has tanh(M) M and cosh(M) « 1, so that we obtain 
M = njj(47rif )/(A2 + A^) and A(iJ) « Ah=o from Eq. 
(10) and Eq. (9), respectively. Since AH=a ^ {Ar,Aa) 
as T < 0.9Tbkt, and A|j, ~ a?/RK, we obtain (using 
Sir'^KkBT/'^l = 1/A): 
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where we recognize fiux expulsion (n = Ba? /^q k, 0) and 
the Meissner effect (— 47rAf = H/ (1 — rff) in the presence 
of a large demagnetization factor rj 1 — d/Ras expected 
in a thin film[17, 18]. At large field instead cosh(A/) w 
e*^/2 and tanh(M) w 1. We then obtain that M « 



log(2A^ ^/p) from Eq. (9), and using A^ 
from Eq. (10) we get: 

^BKT 
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with A{T) = fi/keT - {K/2)log{4TTK). The linear 
regime (12) survives up to a field iJf that can be de- 
termined by the numerical solution of Eqs. (9)-(10). As 
it is shown in Fig. 2b, is very low (~ 10~® G) but 
finite at Tkt- For this reason, the field-independence of 
M at criticality implied by Eq. (13) is only valid above 
Hf, below which M oc —H, as expected. This low-field 
crossing to a linear behavior is missing in Ref. [11] where 
M is calculated as a function of B. However, at large 
fields where B H the dependence of M{B) on log(i3) 
derived there coincides with Eq. (13), apart from an ad- 
ditional B dependence of M at criticality that cannot be 
checked with the present variational calculation. Finally, 
we notice that at T well below Tbkt an estimate of 
can be obtained analytically by matching the high-field 
and low-field solutions for Af at A w Aa: 



if)log(2AA)-log(8i^5)] 



(14) 

which reduces for T ^ to the standard definition 
of first critical field in a SC film, Hci{T ^ 0) = 
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(<I>o/47rAi?) log(A/2a) +47r/i/<J)o^[18]. Indeed, as we can 
see in Fig. la, at low T the magnetization displays a 
sharp kink at i?ci and increases just above it, as indeed 
expected at the threshold of flux penetration (see also 
np m Fig. lb). However, at higher temperatures such a 
kink in M disappears due to thermal smearing and the 
minimum of M is located at a field higher than Hd . 
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FIG. 2: (Color online) (a) M{H) above and below Tbkt 
from the numerical solution of Eqs. (9)-(10). (b) Solid line: 
the threshold field H^''' as a function of T. The points show 
the analytical estimates (14)-(15), which agree with the nu- 
merical result except in a small range near Tbkt- (c) M{H) 
in logarithmic scale (curves are spaced by IK). 



At T > Tbkt, i-e. K < 2, M shows again a crossover 
from a linear to non-linear behavior at a field Hf^. To 
estimate fff we can expand the hyperbolic functions in 
Eq.s (9)-(10) around M = 0. For T sufiiciently above 
Tbkt so that Aj,^^ ~ (4ifg)2/(2-^) » {Aa,Ar) we 
obtain the approximate solutions M ~ [ATrKnu) / A"^ and 
A2 = A]j^^[l + l/{2- K){Hi/^ay], where ^ is the 
zero-field correlation length defined above. When the 
second term in the square brackets is ^ 1 the deviations 
of with respect to A|^^q are neghgible, so that 
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At T sufficiently close to Tbkt screening effects cut-off 
both A (i.e. ^) and so that the estimate (15) is no 
more valid, Hf attains a finite value and merges with 
the field discussed above, see Fig. 2b. As it was 
known[ll, 19] the functional dependence of the low-field 
magnetization M on the BKT correlation length ^ in Eq. 
(15) is the same as in the GL theory[15]. However, the 
critical region H < Hj^ where such a dependence is valid 
turns out to be remarkably smaller than in the standard 
GL theory[15], because ^ diverges much faster than in 
the GL case as r — > Tbkt- 

In conclusion, we proposed a new theoretical frame- 
work to investigate the KT physics of 2D superconduc- 



tors in a finite magnetic field, as given by the modified 
sine-Gordon model (6) and the definitions (3)-(5) of the 
physical quantities as a function of the applied magnetic 
field H (instead of B). As we showed within a varia- 
tional analysis of the model (6), we obtain a clear de- 
scription of the Meissner phase below Tbkt, and an es- 
timate of the threshold field H^'^ for the appearance of 
non-linear effects. Above Tbkt the shrinking of the lin- 
ear regime with respect to standard GL fiuctuations is 
a typical signature of the faster divergence of ^ within 
the BKT theory. These results can shed new light on the 
physics of vortices in cuprates. Indeed, taking into ac- 
count that in layered superconductors the intrinsic cut-off 
Rj is provided by the interlayer coupling J± instead of 
A, Rj ^ ay/ J/J±, our 2D calculations can be applied to 
these systems in all the (T, H) range where A ^ 1/ Rj 
(so that for example large demagnetization effects are 
not expected in layered systems). Thus, the persistence 
of a non-linear magnetization up to 77 ~ 0.01 T in a wide 
range of temperatures above Tbkt found experimentally 
in Ref. [7] can be a signature of the rapid decreasing of 
as T ^ Tbkt, which does not contradict but even- 
tually support the KT nature of the SC fiuctuations in 
these systems. Moreover, since £, increases as ^ increases, 
the extremely low values of i/f measured in Ref. [7] sug- 
gest a value of ^ larger than fj,xY , in agreement with the 
result of Ref. [13] based on the analysis of the superfiuid 
density, and call for a deeper investigation of the normal 
phase existing in the vortex cores. 



[1] V. L. Berezinsky, Sov. Phys. JETP 34, 610 (1972); J. M. 
Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 (1973). 
P. Minnhagen, Rev. Mod. Phys. 59, 1001 (1987). 
D. Matthey, et aL, cond-mat/0603079. 
R. W. Crane, et aL, Phys. Rev. B 75, 094506 (2007). 
A. Pourret, et aL, (2007), cond-mat/0701376. 
Z. Hadzibabic, et aL, Nature 441, 1118 (2006). 
L. Li, et aL, Europhys. Lett. 72, 451 (2005). 
M. A. Cazalilla, A. lucci, and T. Giamarchi, Phys. Rev. 
A 75, 051603(R) (2007). 

S. Doniach and B. A. Huberman, Phys. Rev. Lett. 42, 
1169 (1979). 

P. Minnhagen, Phys. Rev. B 23, 5745 (1981). 
V. Oganesyan, D. A. Huse and S. L. Sondhi, Phys. Rev. 
B 73, 094503 (2006). 

T. Giamarchi, Quantum Physics in One Dimension (Ox- 
ford University Press, Oxford, 2004). 
L. Benfatto, C. Castellani, and T. Giamarchi, Phys. Rev. 
Lett. 98, 117008 (2007). 

1. Nandori, et aL, J. Phys. Cond. Matt. 19, 236226 
(2007). I. Nandori, et aL, arXiv:0705.0578. 
A. E. Koshelev, Phys. Rev. B 50, 506 (1994). 
J. Jose, et aL, Phys. Rev. B 16, 1217 (1977). 
L. D. Landau and E. M. Lifchitz, Electrodynamics of 
Continuous Media (Pergamon, Oxford, 1984). 

A. L. Fetter and P. C. Hohenberg, Phys. Rev. 159, 330 
(1967). 

B. I. Halperin and D. R. Nelson, J. Low. Temp. Phys. 
36, 599 (1979). 



[2 
[3 
[4 
[5 
[6 
[7 
[8 

[9 

[10 
[11 

[12 

[13 

[14 

[is: 

[16 
[17 

[18 

[19 



